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Abstract 

A low temperature model system consisting of a central spin cou- 
pled to a spin-bath is studied to determine whether interaction among 
bath spins has an effect on central spin dynamics. In the absence of 
intra-environmental coupling, decoherence of the central spin is fast 
and irreversible. Strong intra-environmental interaction results in an 
effective decoupling of the central spin from the bath and suppression 
of decoherence. Weaker intra-environmental coupling reduces but does 
not eliminate decoherence. We believe that similar behaviour will be 
observed in any system with a self- interacting environment. 

PACS numbers: 03.65.Yz, 03.67.Lx, 05.30.-d 

1 Introduction 

Exact theories for the dynamics of open quantum systems require solution 
of the Schrodinger equation (or an equivalent formalism) for the full system 
plus environment degrees of freedom. This is impossible to do analytically 
except when the environment consists of a set of independent harmonic os- 
cillators (e.g., Caldeira-Leggett model fLj, spin-boson model @) or spins |§. 
Neglect of intra-environmental coupling is thus strongly motivated by rea- 
sons of mathematical convenience. The purpose of the present study is to 
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determine whether this approximation is justified physically. Our model con- 
sists of a small (i.e., highly quantum) central spin interacting with a bath of 
self- interacting spins representing the degrees of freedom of the environment. 
We shall sometimes refer to this central spin as the subsystem. 

Intuition suggests that coupling between environment degrees of free- 
dom can significantly affect the properties of the bath and consequently, 
through subsystem-environment coupling, the dynamics of the central spin. 
Mutual interactions of the bath modes allow energy exchange without using 
the subsystem as an intermediary. Consider an initial state of the subsystem- 
environment which is perturbed away from equilibrium. With intra-environ- 
mental coupling the rapid initial flow of energy toward a new equilibrium 
state will largely occur within the environment. Without intra-environmental 
coupling the energy must flow through the subsystem, leaving it strongly 
entangled with the bath and hence badly decohered. [Note that the deco- 
herence effects discussed in this manuscript differ greatly from those which 
affect macroscopic subsystems @].] 

Intra-environmental interactions may also alter more abstract properties 
of the bath. The classical dynamics of a bath of oscillators can change from 
regular to chaotic when non-linear interactions are added M. In the quan- 
tum case, such coupling modifies the statistical properties of the energy levels 
and eigenstates of the bath. Specifically, the energy spectra of quantum sys- 
tems with chaotic classical counterpart exhibit level repulsion, while systems 
with regular dynamics in the classical limit show a clustering of the energy 
levels ||. The structure of the bath eigenstates also changes when the dy- 
namics of the bath undergoes a transition from regular to chaotic. This effect 
can be analysed by considering the form of the Wigner functions of energy 
eigenstates: for chaotic systems the Wigner functions spread more or less uni- 
formly over the energetically available phase space J/J, whereas in the case of 
regular systems they are more lumpy. Moreover, these spectral signatures of 
chaos have dynamical consequences M. It is therefore reasonable to expect 
that dissipation and decoherence will assume different forms for coupled and 
uncoupled baths. 

A qualitative understanding of the effects of bath self-interaction on de- 
coherence, while of obvious theoretical interest, might also have important 
applications. Minimization of decoherence is essential for the development 
and implementation of a number of new technologies such as quantum com- 
puting |J, laser control of chemical reactions [|K]] and molecular electron- 
ics |TT|. For quantum computing some proposed physical platforms, such as 
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laser manipulated cold ions in traps in near vacuum |T2] can very effectly 
minimise decoherence. Whether such platforms can be scaled to the 10 5 ions 
needed to perform useful computations like factorisation of large integers is 
unclear |12[]. Proposed solid state platforms such as single-electron quan- 



tum dots embedded in a semiconductor |L3| are readily scaleable HI , but 
decoherence is a serious obstacle. In principle decoherence and dissipation 
can be reduced through judicious choice of states for implementation of the 
qubit fl5 |. (A trivial example would be to avoid states with electric-dipole 
allowed transitions.) Since the solid state in principle provides enormous free- 
dom over the choice of qubit and matrix, further theoretical insight into the 
mechanisms of decoherence might prove extremely useful in selecting optimal 
configurations. 

The effect of environmental self-interaction is almost certainly of impor- 
tance in the solid state. Unfortunately, as we noted above, the analysis of the 
mechanisms of decoherence and dissipation in self-interacting environments 
is a problem that defies exact analytical treatments. Standard approaches, 
like the Feynman- Vernon influence functional method ]16[ or the Nakajima- 
Zwanzig projection technique fl7j cannot be successfully applied unless re- 
strictive hypotheses are made on the nature of the bath and its coupling to 
the subsystem of interest. Thus, for exact results we must rely on what we 
can calculate numerically. Environments with few degrees of freedom should 
in principle provide much of the qualitative information we seek. Such nu- 
merical studies also serve a second purpose, namely to provide exact results 
which can be used as a touchstone to test the reliability of the approximate 
analytical methods which describe open systems (e.g., Redfield SRA [19 
master equations and other approaches p 



In this paper we study numerically and without approximation the dy- 
namics of a spin-1/2 subsystem coupled to a bath of N interacting spin-1/2 
modes. While spin-spin-bath models have been previously studied |3|, |2l| 
intra-environmental coupling has been neglected. Many different physical 
systems can be represented with such a model. Examples include a spin 
chain, a set of magnetic impurities in a solid |TB[ or an array of magnetic mo- 
ments in a molecular crystal. Alternatively, the system can be interpreted as 
a model of an atomic impurity in a crystalline solid at low temperature. The 
impurity is vibronically coupled to a bath consisting of the phonon modes of 
the solid At low temperatures, the probability that the impurity is in 
an electronic state other than the ground state is negligible. However, the 
impurity can be prepared in a chosen (non-radiative) excited state and so is 
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effectively a two-level system. By the same token, in the low-temperature 
limit the phonon modes reduce to oscillators with only two accessible levels 
and can therefore be described as spin-1/2 modes (i.e., a^a — > a z , a^+a — ► a x , 
etc.). Jahn- Teller interactions couple the central spin to the bath. Phonon- 
phonon interactions inter-couple the bath spins. In our study we consider 
both ferromagnetic and antiferromagnetic interactions. 

We show that environmental self-interaction has an extremely important 
effect on decoherence. Strong intra-environmental interactions effectively 
decouple the subsystem from the environment, and even weak interactions 
qualitatively change the manner in which phase information is lost. We see 
no reason to suppose the effect is peculiar to the spin-spin-bath model, and it 
should prove of importance in areas of physics where preservation of quantum 
interference is paramount. 

Section ^explains our model in detail. In section ^ we outline the numer- 
ical methods we employ to solve the subsystem dynamics. Numerical results 
for subsystem dynamics are discussed in detail in section [| for a wide range 
of intra-environmental couplings. We summarise our findings in section ^ 
and discuss their relevance for more general environments. 

2 The spin— spin-bath model 

In mathematical terms, our model is defined by the Hamiltonian 

H = H + H B + H I (1) 

where 

#o = ^?> (2) 
is the Hamiltonian of the central spin (denoted by the superscript 0), 

JV i.. N N-l N 

Hb = Y, ~Y°f + P E + a E E (3) 

i=l i=l i=l j=i+l 

is the Hamiltonian of the bath, and 

N 

#/ = AoE^M 0) (4) 
i=i 
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describes the interaction between the central spin and its environment. For 
simplicity, in the rest of this paper we adopt a system of units such that 
h — 1. Note that the terms proportional to (3 and A in Hamiltonians (|3]) 
represent anharmonic phonon-phonon interactions or non-linear spin inter- 
actions, depending on which interpretation of the model one has in mind. 

To complete the definition of the model, we have to specify the values 
of the various parameters which appear in the Hamiltonian. We assume 
that the frequencies of the bath spins are positive random variables with the 
Debye probability density 



appropriate for the low energy acoustic modes of a crystal. This is equivalent 
to saying that the bath has an Ohmic spectral density with a sharp cut- 
off at uj = uj c which is not necessarily the Debye frequency itself. At the 
temperatures we consider, modes with large frequency are unlikely to be 
populated and so it makes sense to choose a smaller uj c to reflect this fact. 
For the frequency of the central spin, we arbitrarily chose the value ujq = 
0.8288cj c . Note that the results obtained in this manuscript do not depend 
crucially on the specific form (|5|) of the frequency distribution. In fact, we 
repeated our calculations with the box distribution 



and found the qualitative behaviour of the model unaltered by the change. 

To simplify the form of mathematical expressions, we set u c — 1, Ao = 1, 
and we varied the relative strength of the intra-environmental interactions by 
letting the parameter A range in the interval from A = (bath without inter- 
nal spin-spin coupling) to A = 10 (strong bath self-interaction). In addition 
to considering positive values of A, which correspond to antiferromagnetic 
interactions, we investigated the case of ferromagnetic couplings, letting A 
assume negative values in the interval [—10,0]. 

We set (3 = 0.01. An important consequence of the fact that (3, although 
small, is not zero, is that the Hamiltonian (|I]) cannot be reduced to block 
form. To understand this point, we observe that the Hilbert space of the 




for < uj < uj c 
for uj c < uj 



(5) 




for < uj < uj c 
for uj c < uj 



5 



system ([!]) is spanned by the vectors 



|0> 


= 10)^10)^-1.. 


- |0>i[0> 


11) 


= |0)jv|0);v-i.- 


■|0)i|l> 


|2> 


= |0)jv|0)iv-i.- 


■|l)i|0> 


\2 N+1 - 1) 


= |l)jv|l)jv-i • • 


■ ll>lll> 



where the symbols |1), and |0)j denote the 'up' and 'down' states of the i— th 
spin, i.e., the eigenstates of the Pauli z-spin matrix uW with eigenvalues 
+ 1 and —1 respectively. Note that the basis states @ can be conveniently 
interpreted as binary representations of integer numbers ranging from to 

2 N+1 the yector 

\k N , fcjv-i, ■ ■ ■ ,ko) = \kN)N\kN-i)N-i ■ ■ ■ \ko)o (7) 
(with ki — 0,1 for i = 0, . . . , N) is associated to the integer 

N 

k = k + 2k x + . . . + 2 N k N = 2 n fcn- (8) 

n=0 

For (3 = the Hamiltonian ([!]) can be reduced to block form by regrouping 
the basis vectors @ in two sets defined by the condition that the states of 
each set have an even or odd number of spins 'up'. In fact, for (3 = all 
terms of the Hamiltonian (|l]) have the effect of flipping either zero or two 
spins at once, thereby leaving invariant the subspaces spanned by the 'even' 
and 'odd' basis states. The terms proportional to j3, on the other hand, flip 
just one spin and therefore connect the 'even' and 'odd' subspaces, making 
the Hamiltonian ([]]) irreducible. 



3 Numerical approach 

Assume that the bath is initially in thermal equilibrium at temperature T 
and that the central spin is in the excited state 1 1 ) o - [In the case of an 
impurity in a insulating solid such an initial state could be prepared using 
a fast laser pulse with a frequency matching a transition of the impurity 
but lying in the crystal's band gap.] The initial density matrix of the whole 
system therefore has the product form 

p(0) =Po(0)®Pb(0), (9) 
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with 



p Q (0) = |1) 00 (1| and p B (0) = (1/Q) exp(-H B /kT) 
where Q = Ty b [exp(— H B /kT)\. After introducing the notation 

H B \4>W) = E n \4><&) (10) 
for the 2 Ar eigenvalues and eigenvectors of the bath Hamiltonian and 

|^(O)) = |l) o ®|0f ) ) (11) 

for the corresponding initial conditions of the total system, we can write the 
initial density matrix ([|) in the form 

2^ p -E n /kT 

p(o) = £|«)))— — (^n(0)|. 

n=l t> 

To study the dynamics we have determined the states 

|Vn(*)> = exp (-iHt) |Vn(0)) (12) 

which evolve from (|Tl|). Once states ( |12"D are known the evolved density is 
constructed via 

2 N e -E n /kT 

p(t) = exp(-iHt)p(0)exp(iHt) = ]T |^ n (t))— — - (^(t)|, (13) 

n=l W 

and the reduced density of interest is 

-E n /kT 

Po(t) = Tr B [p{t)\ = T,(^ ) \Mt))— F5 —(Mt)\^ ) ). (14) 

We used two complementary techniques to compute the dynamics (|i"2"|). 
For small baths, iV < 11, we numerically diagonalised both the bath Hamil- 
tonian (^|) and the total Hamiltonian ([!]) using standard Householder routines 
(see Ref. P3[). This method gives the exact reduced density at all times and 
allows us to consider a range of temperature. However, it cannot be used for 
baths composed of a large number N of spins since matrices of size 2 N x 2 N 
quickly exceed computer resources. For large baths with iV > 12 spins, we 
exploited the low-temperature limit. For kT — > 0, high-energy eigenstates of 
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the bath are depleted and one can truncate the sum in Eq. flTB| ) to the first 
M terms, with M <C 2^, so that the density matrix takes the form 

M p -E n /kT 

p{t)^Y,\Ut))—F—(Ut)\- (is) 

n=l W 

Thus a complete diagonalisation of the bath Hamiltonian becomes unnec- 
essary. We chose a bath temperature of kT = 0.02 for which the number 
of terms M needed in (O) is 20. We checked that neglected terms were 
irrelevant by evaluating the probability that the bath be in an eigenstate of 
energy E > E^o- For baths with N < 11 spins with < A < 10 we obtained 
an upper bound 

20 

P(E B > E 20 ) = 1 - E e~ EnlkT IQ < 10- 4 , 

n=l 

showing that bath eigenstates with E > E20 are unpopulated. As a further 
check, for baths composed of more than 11 spins we estimated the ratio 
R = p 20 / Pl = exp( J B 1 - E 20 ). For iV = 14, < A < 10 and kT = 0.02 the 
probability ratio was less than 10~ 6 . 

To determine the M bath eigenstates of lowest energy, we used ARPACK 
routines based on the Lanczos algorithm(see [f24[). Evolved states ([i~2|) were 



calculated using a Runge-Kutta algorithm of eighth order [25]. Neither the 
programs for the partial diagonalisation of the bath Hamiltonian, nor the 
Runge-Kutta subroutine required that the whole Hamiltonian matrix be 
stored in the computer memory, but only that the matrix-vector product 
H\ip) be defined. This more efficient use of computer resources allowed us to 
consider baths of up to 14 spins. 

We calculated H\ip) given an input state \ip) by iterated calls to subrou- 
tines which multiplied by and Consider multiplication by g%> as an 
example. With vectors (H) as a basis for our Hilbert space 



a 



(')i 



2 jv+i_ 1 



E (16) 



k=0 



thus reducing our problem to that of finding an efficient way to multiply the 
basis vectors (|]) by erg'. Since the matrix a^> has the effect of flipping the 
i-th spin, one has 



a 



{ ^\k N , . . . , hi, . . . , ho) = \k N , . . . , hi, . . . , k ) 
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where ki = 1 if ki = and fcj = if ki = 1. Multiplication by therefore, 
replaces the fe-th component of \ip) (where k = k + . . . + 2 l ki + . . . + 2 N k^) 
with the fc'-th component (where k' — ko + • • • + 2*fcj + . . . + 2 N k^) and vice- 
versa. In binary representation, the numbers k and k' differ by a single bit 
(the i-th bit) and one can therefore obtain k! from k using Fortran intrinsic 
functions. Specifically, we used the XOR-function (exclusive or) to flip the 
i-th bit of the k-th. state. Multiplication by a® can be similarly implemented. 

After determining the evolved density (0), we traced out the bath degrees 
of freedom to obtain the reduced density ([14]). As indicators of quantum 
coherence, we chose the polarisation and entropy of the central spin defined 
respectively as 

P(t) = Tr[p (t)a] (17) 

and 

1 /i_p2\ p n + p\ 
So(t) = -Tr [p (t) \np (t)} = -- In [^^J " 3 ln [t^p) ' (18) 

Here P = \P\ denotes the modulus of the polarisation vector P. Note 
that (|17D contains as much information as the reduced density itself. In 
fact ( |T4| ) can be expressed in terms of P via 

p (t) = 1 -(l + P(t).a). 

4 Dynamics of the central spin 

Here we examine the effects of intra-environmental couplings on the dynamics 
of the central spin. We consider the antiferromagnetic and the ferromagnetic 
cases separately. High temperature results apply only to the case where the 
spin-bath represents true physical spins. 

4.1 Antiferromagnetic interactions 

To evaluate the effect of antiferromagnetic interactions we calculated the 
dynamics of the central spin for values of A ranging from A = (uncoupled 
spins) to A = 10 (strong coupling). As a point of reference, note that in the 
absence of subsystem- environment coupling P z {t) - initially one - undergoes 
periodic fluctuations to slightly smaller values. The components P x {t) and 
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P y (t) - initially zero - oscillate about zero with the same period and similar 
small amplitude. 

The entropy fllBf ) is shown in Fig. [I]. Figs. |2], and || show components 
of the polarisation (|T7|). In the figures, each curve corresponds to a 
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Figure 1: So(t) versus u c t (14 bath spins). 



different value of the intra-environmental coupling A (reported in parenthesis 
to the right). The thermal average was computed for a bath at temperature 
kT = 0.02. The results shown are for N = 14 but are also representative of 
smaller baths with even numbers of spins (i.e., 8, 10 and 12). The time scale 
considered is long enough for the system to attain its asymptotic condition 
as can be seen from Fig. |5] which shows the entropy of the central spin on a 
much longer time scale (data obtained for iV = 10 spins). 

When A = (i.e., bath spins non- interacting), the central spin under- 
goes rapid decoherence, with the polarisation quickly falling to zero and 
the entropy So(t) simultaneously approaching the maximum value S™ ax = 
ln(2) ~ 0.693147.... As A increases, however, the entropy So(t) tends to 
progressively smaller asymptotic values and the components of the polarisa- 
tion vector approach the dynamics of a central spin evolving in isolation. In 
other words, the existence of a strong interaction among bath spins suppresses 
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Figure 2: P z (t) versus cu c t (14 bath spins). 
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Figure 4: P y (t) versus cu c t (14 bath spins). 




Figure 5: So(t) versus uj c t (10 bath spins). 
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environmental decoherence. 

This apparently puzzling result is a straightforward consequence of the 
fact that strong interactions between bath spins produce an almost com- 
plete decoupling of the central spin from the bath. This can be verified by 
considering the thermal average of the interaction Hamiltonian (§), defined 
as 

2 N 

(Hj(t)) = J2(Mt)\Hi\Mt))e~ E " /kT /Q. (19) 

n=l 

Note that (Hi(t)), rather than Ao, is the physically relevant quantity deter- 
mining the strength of the interaction. In fact, the interaction may be small 
even if A is large. The evolution of (Hi(t)) is displayed in Fig. |6|, which 
shows that as A increases the effective interaction of the central spin with 
the bath tends to zero (data obtained for 10 spins). Quantitatively, (Hj) 
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Figure 6: Thermal average (Hj) of the interaction Hamiltonian versus u c t 
(10 bath spins). 



averaged over the time interval [0 : 2500u; c x ] decreases from (Hi) = —3.07 
for A = to (i/j) = -0.01 for A = 10. 

There is a simple explanation for the observed behaviour of the average 
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interaction term fllPp. Define the total bath spin 



N 

£ = £a«, (20) 
i=i 

(i.e., E x = Y,iLi cr x' > etc.) and rewrite the interaction Hamiltonian (fD and 
bath Hamiltonian (|3[) in the suggestive forms 

#, = Aq^E, 

and 



2 



A' 



X 



+ /3E, + E>0O. (21) 
i=i z 

Eq. (|2l|) shows that for A ^> c<j c the bath Hamiltonian takes the approximate 
form 

A r 



H B = - 
2 



Nl 



+ 0(A U ) (22) 



which is essentially proportional to Y? x . Thus for large values of A the bath 
eigenstates ([10]) must be approximate eigenvectors of T, x and those of lowest 
energy must correspond to the eigenstates of E x with zero eigenvalue. (That 
many such eigenstates exist is a consequence of the fact for every value of the 
macroscopic variable E^ there are many corresponding microstates). Thus, 
for low temperature and large A the relevant diagonal and off-diagonal matrix 
elements of the environment coupling operator S x will be zero. 

These conclusions are confirmed by the numerical evaluation of the expec- 
tation values of on the bath eigenstates fllQ|) . Fig. |7| reports {4>^\H x \(j)^) 



as a function of the index n which orders the eigenstates \(f)^) in ascending 
energy (data obtained for N = 10 so n runs from 1 to 1024). Data are shown 
for A = and A = 10. For A = the expectation values of S^. are small but 
non-zero (this is more evident in the inset, which displays the data for the 
100 lowest energies). For A = 10 the expectation value of T, x has a sort of 
step-like behaviour and attains relatively large values for high-energy bath 
eigenstates. At low temperatures only low-energy states matter and for these 
(^lEal^f)) is orders of magnitude smaller than for the A = case (see 
inset in Fig. 0). This effect can be quantified by computing the quantity 

i 20 

(^) = ^E(4 B) |S,|0i B) ). (23) 

ZU n=l 
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Figure 7: Expectation value {<^'\ll x \<^') versus n for A = and A = 
The stacked inset shows on a larger scale the first 100 values. The data 
obtained for a 10-spin bath. 
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One finds 



(E»> 



-0.161951 for A = 
-6.0 • 10" 5 for A = 10 



(24) 



for N = 10 spins which is indeed small for strong coupling. 



4.1.1 High temperature limit 

Fig. [7], for A = 10, shows that the expectation values of S x are large for high- 
energy bath states. Populations in these states are zero at low temperatures 
but increase with temperature. The thermal average interaction ([B5J) will 
therefore also increase with temperature, effectively coupling the central spin 
to its bath. As a consequence, one can expect that the self-interacting bath 
will behave more and more like an ordinary bath of uncoupled spins when 
the temperature is raised. Fig. H shows the time behaviour of the central 
spin entropy for a bath of 10 coupled spins (with A = 10) at various tem- 
peratures. Entropy increases monotonically with temperature, progressively 
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Figure 8: S versus u> c t for A = 10 for various temperatures of the bath (10 
bath spins). 



approaching the behaviour characteristic of the uncoupled bath. However, 
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even at high temperature (kT = 300) the entropy remains lower than for 
a bath of uncoupled spins: this shows that even at high temperatures the 
spin-spin coupling has a reductive effect on decoherence. These features are 
confirmed in Fig. |9|. Raising the temperature lowers the asymptotic value of 
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Figure 9: P 2 versus uj c t for A = 10 for various temperatures of the bath (10 
bath spins). 



P z (t) which however never becomes zero. 

In summary, at low temperatures the strong intra-environmental inter- 
actions force the bath spins to align in an antiferromagnetic state of zero 
moment and behave like a single giant spin decoupled from the central spin. 
As the temperature is raised thermal fluctuations eliminate the spin align- 
ment and switch on the subsystem-environment coupling, making the self- 
interacting bath behave more like an ordinary bath of uncoupled spins. 

4.1.2 Odd numbers of spins 

So far we have considered the behaviour of a central spin coupled to a bath 
composed of an even number of spins. To complete the discussion of antiferro- 
magnetic interactions, we briefly mention a last feature of the self-interacting 
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spin bath 
of spins. 
Figs 



namely modifications which appear for baths with an odd number 



TP] and [XT] show So(t) and P z (t) for N = 11. Similar behaviour 

and 9 spins. Only A = 0, 1 and 10 are 
1 overlap. The data show 



was obtained with baths of 5, 7 
shown because the curves for A greater than 




Figure 10: Sq versus u c t (bath composed of 11 spins). 



that, regardless of whether the number of bath spins is even or odd, S (t) 
is a decreasing function of the intra-environmental coupling. However, for 
iV odd, the average value of the entropy decreases less than for N even and 
oscillates more (compare Figs. [TO] and [I]). Polarisation behaves similarly. As 
shown in Fig. [11], when the intra-environmental interactions are strong, P z (t) 
oscillates around the value P z (t) — 0.8 for baths with odd number of spins, 
whereas one has P z {t) — 1 for a bath with an even number of spins. Nor 
do the oscillations die on a longer time scale: in fact, numerical calculations 
performed via the exact diagonalisation of the Hamiltonian show oscillatory 
behaviour persisting on time scales an order of magnitude longer than that 
in Figs. [Tl] and O. Thus intra-environmental interactions have the effect 
of suppressing decoherence regardless of how many spins compose the bath; 
however, this effect is more marked if the bath is composed of an even, rather 
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Figure 11: P z versus u c t (bath composed of 11 spins). 



than an odd, number of spins. 

The oscillations of entropy and polarisation are also displayed in the 
strength of the interaction between the central spin and its bath. Fig. [1^ 
shows the thermal average of the interaction Hamiltonian @. Comparing 
Fig. [12] with Fig. ffl, it is evident that the interaction between the central spin 



and bath is weaker for N even than N odd. Computing the parameter (|23|) 
for a bath of 11 spins gives 



-0.19180 for 
-4.97 ■ 10~ 3 for 



A 

A = 



= 
10 



(25) 



Comparing Eq. ( p4|) with Eq. (|25|) shows that (E x ) is a decreasing function 
of A for baths of both even and odd numbers of spins but (E x ) for A = 10 
is two orders of magnitude larger for the 11-spin bath than for the 10-spin 
bath. Table [I] shows that as N increases (E x ) (for A = 10) decreases. The 
even-odd difference is thus a finite-size effect which should vanish in the 
thermodynamic limit. 



19 




50 100 150 200 



Figure 12: Thermal average (Hi) of the interaction Hamiltonian versus u c t 
(baths of 11 spins). 



Number of bath spins 


(E x ) for A = 10 


7 


-0.011147 


9 


-0.006724 


11 


-0.004967 



Table 1: (S^) | a=io as a function of N 
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4.2 Ferromagnetic interactions 



The entropy for the ferromagnetic case is shown in Fig. [13|. The data were 
obtained for a bath of 10 spins; we use a semi-logarithmic scale to distinguish 
the curve for A = —2 from the x-axis. We do not show data for values of 
A < —2 because they overlap. The behaviour of the entropy is similar in the 
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0.001 



1e-05 



1e-06 




Figure 13: In (So) versus u c t (10 spin bath). 



ferromagnetic and antiferromagnetic cases; the only significant difference is 
that in the ferromagnetic case a suppression of the entropy can be achieved 
with weaker intra-environmental interactions. 

The most relevant difference between antiferromagnetic and ferromag- 
netic baths emerges when one considers the behaviour of the polarisation 



vector (17) shown in Figs. 14, 15, and Gil We have chosen a time scale 



appropriate to the fast oscillations of the polarisation components (but long 
enough to be representative of the long-time behaviour of P{t)). For A < —2 
the components of the polarisation vector assume a strongly oscillatory be- 
haviour. This is consistent with the suppression of entropy since entropy 
depends only on the norm of the polarisation vector which can stay close to 
unity even if individual components of P oscillate in time. Such polarisation 
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Figure 14: P z versus cu c t (10 spin bath). 
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Figure 16: P y versus u c t (10 spin bath). 



dynamics, however, is very different from that observed in the antiferromag- 
netic case where P z {t) — 1 and P x (t) ~ P y (t) ~ at all times. Term (19) 
is not negligible for ferromagnetic interactions, as can be seen in Fig. [17], 
where we represent the thermal average of the interaction Hamiltonian as a 
function of time. For A < —2, the interaction term (Hj(t)) becomes almost 
independent of the strength of the spin-spin coupling and oscillates around 
a non-zero value. 

This behaviour can be explained by writing the bath Hamiltonian in 
the form (|2T[) and by observing that in the limit of strong ferromagnetic 
interaction 

|A| 



H B = -— 
2 



s: - ni 



0(A°). (26) 



Note that the antiferromagnetic and ferromagnetic Hamiltonians (|22| ) and ( |26[ 
are identical but with opposite sign. The bath eigenstates of lowest energy 
are eigenstates of E x in both cases, but with (<p^\H x \<p^) ~ in the anti- 
ferromagnetic case and \((p^\J2 x \(f)^)\ ^> 1 in the ferromagnetic case. This 
is confirmed by Fig. [18], which represents ((j)^\T, x \cj)^) as a function of the 
index n. Not surprisingly, the distribution of points for A = —10 is almost 
the reverse of the distribution for A = 10 shown in the corresponding Fig. [?]. 
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Consequently the coupling term ( |T§| ) is significantly larger for a ferromagnetic 
bath than it is for a antiferromagnetic one. Thus, diagonal matrix elements 
of the coupling operator S x are large while off-diagonal matrix elements are 
again zero. 

It is important to note that the oscillations of the average interaction 



Hamiltonian (19) are synchronous with those of P x (t): more precisely 

P'(t) = (H x (t))/P m (t) 
is practically constant in time and, for A < —2, it takes the value 

P' ~ -10 

regardless of the strength of the spin-spin coupling. The value of /3' is eas- 
ily explained. The matrix element (<j)f\Y} x \(j)f) corresponding to the lowest 
energy bath eigenstate is -10 while that of the next lowest is 10. How- 
ever, the energy gap between these states is .2 leading to a population ratio 
Vilvi ~ 10~ 5 at kT = .02. Hence, the lowest energy state is the only one 
populated (in sharp contrast to the antiferromagnetic case) and its S x eigen- 
value is -10. This also explains why decoherence is more easily suppressed in 
the ferromagnetic case. 

Thus the central spin interacts with the bath through 

tf 7 e// = /3Vf, (27) 

and the evolution of the central spin is dictated by Hamiltonian (0) with a 
renormalised (3 parameter: 

-+ p = p + ~ -10. (28) 

Solving the Heisenberg equations for the renormalised central spin with initial 
conditions P z (0) = 1, P x (0) = P y (0) = gives 

PS) = ^(1-cosnt) 

P y (t) = sin nt (29) 

p z {t) = l--^-(i-cosm) 
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where f2 = yujQ + A/3 2 . Comparing Eq. fl29|) with the actual data shows 

that ( |2"T| ) is correct. Figs. and |TB| do not show P(t) predicted by 

Eq. (^) because of overlapping of the various curves. 

We therefore conclude that the only effect that a low temperature bath 
with internal ferromagnetic interactions has on the central spin is a renormal- 
isation (^) of the /3 parameter in the Hamiltonian (§) (i.e., a Lamb shift). 
In both the antiferromagnetic and the ferromagnetic cases, therefore, when 
intra-environmental interactions are strong the dynamics of the central spin 
is almost autonomous from that of the bath and is dictated by a Hamiltonian 
of the form 

As a last remark on the ferromagnetic case, we observe that the dynamics 
of the central spin is unaffected by whether the bath is composed of an even 
or odd number of spins. 



5 Summary and Conclusions 

In this work we considered a spin 1/2 subsystem coupled to a low-temperature 
bath of interacting spin-1/2 modes. We focussed attention on the effects of 
antiferromagnetic and ferromagnetic intra-environmental interactions on de- 
coherence. In both cases strong intra-environmental interactions suppress 
decoherence by making the dynamics of the central spin almost autonomous 
from the bath itself. More precisely, strong antiferromagnetic couplings 
among bath spins make the average value of the subsystem-environment in- 
teraction Hamiltonian vanish, thereby making the central spin evolve with 
its unperturbed and uncoupled Hamiltonian. Strong ferromagnetic couplings 
among bath spins, on the other hand, cause a Lamb shift of the subsystem 
Hamiltonian but otherwise leave it to evolve in isolation. 

These effects can be schematically explained by considering with Hj = 
SB where S and B are subsystem and bath operators. Now, the eigenstates 
of Hb are also eigenstates of B for strong intra-environmental coupling in 
our model. Hence the off-diagonal matrix elements of B in this eigenbasis 
are zero and so @) cannot couple bath eigenstates. Since the initial states of 



the full system are proportional to bath eigenstates for initial conditions ( |TT| ) 
it follows that the time evolved reduced density must be of the form 

p -E n /kT 

pQ (t) = £ e -i(Ho+SB n )t pQ ^y(H 0+ SB n )t (3Q) 

n W 
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where B n is the eigenvalue of B corresponding to eigenstate \4>^) of H B . In 
the antiferromagnetic case we found that all relevant B n were zero for strong 
intra-environmental coupling and hence the dynamics was free of decoher- 
ence. For the ferromagnetic case we found that only the lowest energy state 



was populated at low temperature and so only one term contributes to (pCf ) 
and again the dynamics is coherent, but with a Lamb shift. 

We expect similar effects to occur in more general baths of coupled an- 
harmonic oscillators. The key issue as we have seen is what happens to the 
matrix elements of the spin-bath coupling Hamiltonian, in the eigenbasis of 
the bath, when bath self-interactions are turned on. To suppress decoher- 
ence the off-diagonal matrix elements must be small. Clearly in an integrable 
bath (i.e., no self-interaction) some off-diagonal matrix elements will be large 
due to selection rules. One would thus expect strong decoherence for inte- 
grable environments. In the case of chaotic environments it is known that 
the off-diagonal matrix elements are of order h N ~ x smaller than the diago- 



nal matrix elements [ 26fl , where N is the number of environmental modes. 
Since Planck's constant is small and N is very large, the off-diagonal matrix 
elements for a chaotic environment are negligible. Hence, we expect at most 
a Lamb shift of the subsystem for each eigenstate of the bath and hence a 
reduced density like (|30D . If in addition the diagonal matrix elements (i.e. 
B n ) vary slowly with energy then at low temperatures we should obtain a 
reduced density 

p (t) = e -i(H +SB )t po ( y(H 0+ SBo)t (31) 

where B Q is a representative low energy diagonal matrix element of the cou- 
pling operator. Dynamics should therefore be coherent but (possibly) Lamb- 
shifted for systems interacting with low temperature chaotic environments. 

Finally, returning to the specific case of an atomic impurity in a crystal at 
low temperature, consider the limit of strong phonon-phonon coupling. Since 
the Wigner functions of the bath eigenstates are nearly uniform over the en- 
ergetically available phase space J7|, energy is distributed over an enormous 
number of phonon modes. At low temperature, displacements from equi- 
librium of any phonon mode must therefore be small. Since coupling of an 
impurity to a phonon is through its displacement coordinate [p2 |, this cou- 



pling will also be small. Thus, we expect diagonal and off-diagonal matrix 
elements of the Jahn- Teller interaction to be small and decoherence to be 
minimal. Impurity- crystal configurations which might have these attributes 
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include noble gas substitutional or interstitial impurities in diamond (which 
has strong phonon-phonon interactions [f2~7|j). 

The authors gratefully acknowledge the financial support of the Natural 
Sciences and Engineering Research Council of Canada. 
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